arXiv: 1501.02867vl [astro-ph.SR] 13 Jan 2015 


Accepted for publication in The Astrophysical Journal 

Preprint typeset using I4T^ style emulateapj v. 5/2/11 


MODELING MULTI-WAVELENGTH PULSE PROLILES OP MILLISECOND PULSAR PSR B1821-24 

Yuanjie Du‘, Guojun Qiao^, Ping Shuai', Xiaomin Bei', Shaolong Chen', Linzhong Fu', Liangwei Huang', Qingqing 
Lin', Jing Meng', Yaoiun Wu', Hengbin Zhang' Qian Zhang', and Xinyuan Zhang' 

Accepted for publication in The Astrophysical Journal 

ABSTRACT 

PSR B1821—24 is a solitary millisecond pulsar (MSP) which radiates multi-wavelength pulsed photons. It 
has complex radio, X-ray and 7-ray pulse profiles with distinct peak phase-separations that challenge the tra¬ 
ditional caustic emission models. Using the single-pole annular gap model with suitable magnetic inclination 
angle (a = 40°) and viewing angle (C = 75°), we managed to reproduce its pulse profiles of three wave¬ 
bands. It is found that the middle radio peak is originated from the core gap region at high altitudes, and the 
other two radio peaks are originated from the annular gap region at relatively low altitudes. Two peaks of 
both X-ray and 7-ray wavebands are fundamentally originated from annular gap region, while the 7-ray emis¬ 
sion generated from the core gap region contributes somewhat to the first 7-ray peak. Precisely reproducing 
the multi-wavelength pulse profiles of PSR B1821—24 enables us to understand emission regions of distinct 
wavebands and justify pulsar emission models. 

Subject headings: pulsars: general — pulsars: individual (PSR B1821—24) — radiation mechanisms: non- 
thermal — acceleration of particles 


1. INTRODUCTION 

PSR B1821 -24 (PSR B1821-24A or PSR J1824-2452A) 
is a solitary energetic millisecond pulsar (MSP) firs t discov¬ 
ered in the globular cluster M28 (iLvne et ^I1987^ . It has 
a spin period P = 3.05 ms and a period derivative P = 
1.62 X 10“^® ss“^ accor ding to the ATNL Pulsar Cataloguj^ 
(catalogue version: 1.50: [Manchester et al.ll2005h . Due to its 
largest intrinsic P (lPhinnevlll993[) . PSR B1821—24 is be¬ 
lieved to be the most energetic MSP up to date. 

PSR B1821—24 has detectable multi-wavelength emission 
from radio to 7-ray band. It was the first MSP with X-ray 
pulsed emission clearly detec ted by Advanced Satellite for 
Cosmology and Astrophysics (iSaito et al.lll997l) . The X-ray 
pulsed emission of PSR B1821—24 has also been observed 
by other X-ray telescopes. At the radio band, PSR B1821—24 
was the first MSP which ha d been observed to go thr ough a 
microglitch in 2001 March (ICognard & Backe^l2004ft . This 
could be a useful hint to the conventional glitch theory which 
should work for both young and old millisecond pulsars. Sim¬ 
ilar to other MSPs (e.g., PSR J0218-b4232, PSR B1937-b21 
or B1957-1-20), PSR B1821—24 has radio giant pulses with 
a steep spectral index and a high degree of polarization, the 
phase positio ns of these giant pu lses are overla pped with its 
X-ray pulses (iKnight et al.ll2006h . In addition. iBilous et alJ 
(1201^ first presented a 720-2400 MHz phase-resolved radio 
measurements including polarization angle curve, spectral in- 
di ces and giant pulses. 

iJohnson et akl (120131) showed that 7-ray, X-ray and ra¬ 
dio pulse profiles with good signal-to-noise ratio for 
PSR B1821—24, which challenge the pu lsar caustic models 
(iCheng et al.lllOSdiDvks & Rud^l2003h . Such models can¬ 
not self-consistently explain the complexity of its pulse pro¬ 
files. PSR B1821—24 was detected to have strongly linearly 
polarized radio emission. With the coherently dedispersed 
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polarization data, iStairs et aP (119991) used the rotating vec¬ 
tor model (RVM) to fit its polarization curve, and the incli¬ 
nation angle a = 40? 7 ± 1?7 as well as the impact angle 
(3 = 40° ± 10° are derived, note that the viewing angle is 

= a + p. 

Up to date, the multi-wavelength pulse profiles of 
PSR B1821—24 have not been well simulated by cur- 
rent emission mod els due to its unprecedent complexity. 
iVenter et al.l (120121) speculated that PSR B1821—24 have 
plausible phase-aligned radio. X-ray and 7-ray pulse profiles 
which had not been precisely simulated. Lurthermore, this 
pulsar does not have all peaks occurring at th e same phase, 
accord ing to the observed results presented by IJohnson et alJ 

dMl. 

In this paper, the three-dimension annular gap model is used 
to study multi-wavelength (7-ray, X-ray and radio) pulse pro¬ 
files for PSR B1821—24. In Section 2, we will briefly intro¬ 
duce the definition of annular gap and core gap; the multi¬ 
wavelength pulse profiles for PSR B1821—24 will be sim¬ 
ulated, and the altitudes of multi-wavelength emission will 
also be identified. The discussion and conclusion are shown 
in Section 3 and Section 4, respectively. 

2. MODELING MULTI-WAVELENGTH PULSE 
PROLILES LOR PSR B1821 -24 IN THE ANNULAR 
GAP MODEL 

2.1. Annular Gap and Core Gap 

The polar cap region of a puls ar can be divid ed into the an¬ 
nular gap and core gap regions (iDu et alJl2010t) . In the vicin¬ 
ity of the light cylinder, charged particles could not co-rotate 
with the magnetosphere and they escaped. The pulsar must 
output the particles to compensate its loss. Therefore, a huge 
acceleration electric field (i?||) with opposite signs are gen¬ 
erated in both annular gap and core gap regions. This is the 
basic physical picture for the generation of i?|| in the annular 
gap model, which is suitable for both millisecond and young 
pulsars with short spin periods. The detailed calculations on 
i?ll were presented in (iDu et alJLOfll 120131) . A number of 
primary particles are accelerated to relativistic energies by Eu 
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Fig. 1.— Modeled 7 -ray, X-ray and radio pulse profiles for PSR B1821—24. Left panels: multi-wavelength observed pulse profile data, radio and X-ray data 
taken from figure 5 of IJohnson et all t2013l) . Middle panels: two-pole photon sky-map for an inclination angle of a = 40°. Right panels: the corresponding 
simulated pulse profiles (thin black solid lines) cut by only one pole emission with a viewing angle of = 75°. The single-pole annular gap model can well 
reproduce the observed multi-wavelength pulse profiles (thick red solid lines in the right panels). 

(A color version of this figure is available in the online journal.) 


in the annular and core gaps. Then the primary particles em¬ 
anated a great number of 7 -ray photons, most of them are 
converted into positron/electron pa irs by the proce sses of 7 - 
B and photon-photon annihilation (iDu et al.ll2010l) . Because 
of electron/positron pair creation, the plasma in the magneto¬ 
sphere could not be fully charge-separated. 

2.2. Multi-wavelength Observations 

When folding 7 -ray or X-ray pulse profiles of a radio- 
loud pulsar from Termt-LAT data, a good radio timing so¬ 
lution is required by the TEMP02 tim ing software with some 
relevant plugins (iHobbs et al.ll2006h . For PSR B1821—24, 
[Johnson et al.l (120131) presented a precise timing solutiorQ 
which is suitable for generating both 7 -ray and X-ray pulse 
profiles. Using Tempo2 with the “Fermi” plugin or “photons” 
plugin, we were able to reproduce the 7 -ray or X-ray pulse 
profile of PSR B1821—24. To ensure the validity of the ob¬ 
served pulse profiles, we dire ctly took the 7 -rav a nd X-ray 
profile data from Figure 5 of IJohnson et al.l (120131) (see the 
top and middle figure in the left panel of Figure [T]i. The 1.4 
GHz radio profile data was taken from the webpage of “Table 
of Published Ephemerides” of Fermi Science Support Cen¬ 
ter. The observed pulse profiles of these three bands are then 
renormalized and suppressed to compare with the modeled 
results by the annular gap model. 

The radio profile has three peaks (denoted as peak 1, peak 
2 and peak 3 in turn from left to right), while the X-ray and 
7 -ray profiles have two main peaks (denoted as peak 1 and 
peak 2 from left to right) with different peak separations. The 

^ The timing model (par file) of PSR B1821—24 can be found at 
http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/ 


first radio peak is almost aligned in phase with that of X-ray 
and 7 -ray. However, the other two radio peaks have com¬ 
plex phase separation wit h comparison of the s econd peak of 
X-ray and 7 -ray pulses. JJolmson^alj j20 1 31) attempted to 
use the outer-g ap model fC heng et al.lll986l) ^d the two-pole 
caustic model (iDvks & RudakI120031) to explain the complex 
features of multi-wavelength pulse profiles, but failed, then 
they drawn a conclusion that it would be difficult to model the 
actual pulse profiles across all the three wavebands. Instead, 
we will jointly model these three-waveband profiles in our an¬ 
nular gap model, and reproduce their main features including 
peaks and relative phase separation. 

2.3. Modeling Method 

§3.1 of iDu et ^ (I2OI ll) presented the detailed method of 
modeling pulse profiles. The parameterized gap width and 
emissivities were used, which is similar t o the method us ually 
used in the polar cap and slot gap models (lHardingl2013l) . For 
both of the annular gap and core gap, we project all radiation 
from probable emission spots on open filed lines in the emis¬ 
sion region to the resting sky. Meanwhile, both the related 
aberration effect and retardation effect are considered. An as¬ 
sumption of Gaussian emissivities along an open field line is 
adopted in the modeling to speed up calculation, since this 
numerical assumption is consistent with physically calculated 
spectra in the an nular gap model, as presented in Figure 8 of 
iDu et al.l (I2OI ll) . Note that emissivities along different open 
field lines with a varying magnetic azimuthal ( 1 / 1 ) could be 
different, because the out-flowing particles on different open 
field lines (magnetic tubes) with different ip could have dif¬ 
ferent Lorentz factors and particle numbers. This results in a 
patch-like emission pattern (photon sky-map). 












































































































Modeling Pulse Profiles of PSR B1821-24 


3 


TABLE 1 

Model parameters oe multi-wavelength light curves eor PSR B1821 -24 
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The model parameters for modeling radio, X-ray and 7 - 
ray pulse profiles of PSR B1821—24 are listed in Table 1. 
The model parameters are simply introduced below, a and 
( are magnetic inclination angle and viewing angle; k and 
A are two geometry parameters used to scale the peak al¬ 
titude in the annular gap; e is a ratio for the peak altitude 
in the core gap in units of the peak altitude in the annular 
gap; 'ipcut is a special magnetic azimuthal used to constrain 
different regions which Gaussian emissivities with different 
parameters are adopted; ctai and (Ta2 are length scales in 
units of i?LC for the emission region on open field lines of 

— 180° < ips < i’cnt and ^ fps < 180° in the an¬ 
nular gap, respectively; ac is similar as ctai, but a length 
scale for the core gap; erg ai and erg a 2 as well as erg ci 
and CTe_c 2 are the transverse bunch scale for field lines of 

— 180° < ips < fpeut and ip°^^ < ips < 180° in the annu¬ 
lar gap as well as core gap, respectively. T he more detaile d 
description of these symbols is presented in (iDu et al.ll201 ll) . 
We took the values of a = 40° and ( = 75° which are con¬ 
sistent with those of the radio polarization measurement of 
PSR B1821-24 (iStairs etal.|[T9^ . a and C are the most im¬ 
portant parameters in our model. When they are chosen, other 
model parameters which stands for the emission regions are 
then carefully adjusted to reproduce the observed light curve 
of each band. 

2.4. Modeling Results 

The modeled multi-wavelength pulse profiles of 
PSR B 1821—24 are well reproduced by the single-pole 
annular gap model, as presented in Figure 1. The correspond¬ 
ing emission heights for each band are shown in Figure 2. 
From modeling of the radio pulse, it is found that the radio 
emission has a non-caustic origin because of its Gaussian 
emissivities (see Table 1), while caustic emission tends to 
have a constant emissivity along an field line which is a 
conven tional assumption as used in the two-pole caustic 
model (iDvks & Rudakll200^ . The peak 1 and peak 3 of radio 
pulse originate from two regions in the annular gap with 
intermediate altitudes of 0.64i?LC to 0.85i?LC 0.52i?LC 

to 0.78i?LC. respectively; while radio peak 2 is generated 
in the annular gap region with higher altitudes of 0.87i?LC 
to 1.07i?LC- For X-ray pulse, it is generated in the annular 
gap region located at only one magnetic pole. The emission 
heights of X-ray peak 1 are nearly overlapped with those 
of radio peak 1, while the emission region of X-ray peak 2 
is 0.76i?LC to 0.98i?LCj slightly higher than the heights of 
radio peak 3. For 7 -ray pulse, the emission is also generated 
in the annular gap region. Peak 1 is formed by two distinct 
regions with altitudes of 0.36i?LC to 1.20i?Lc; while peak 
2 originates from a region with altitudes of 0.35i?LC to 
0.64i?LC- 

3. DISCUSSION 

An observed integrated pulse profile is just a long-time ac¬ 
cumulated photon counting (for X-rays or 7 -rays) or signal 



Fig. 2.— Emission heights vesus phase for radio, X-ray and 7 -ray pulse 
profles of PSR B1821—24. It is found that shifted emission altitudes for 
each waveband in the single magnetic-pole, which leads to a patch-like 
emission pattern, can explain the complex properties of pulse profiles for 
PSRB1821-24. 

power measurement (for radio emission) based only on a cer¬ 
tain orientation (i.e., line of sight). It is obviously difficult 
to obtain the whole emission pattern on the sky for a pulsar. 
To maximally reproducing the multi-wavelength profiles of 
PSR B1821—24, the values of the most important parame¬ 
ters are chosen as a = 40° and p = 75°, which are consis¬ 
tent with those derived from the radio polariz ation measure¬ 
ment with the RVM. As noted by|CraiJ(l20l3), the effects of 
pulsar magnetic field configuration, emission altitude, orthog¬ 
onal mode jumps, open zone radius and interstellar scatter¬ 
ing would significantly change radio polarization properties, 
which results in the inaccurate values of a and p obtained by 
RVM. Particularly, the situation will be even more complex 
when considering the prop agation effects for a radio photon 
in a pulsar magnetosphere (IWang et al.ll2014f) . As a conven¬ 
tional meth od used by the two-pole caustic model an d outer 
gap model (iDvks & RudakllT003t ICheng et al.ll200()h . a and 
P are are regarded as free parameters, and obtained via fit¬ 
ting pulse profiles. Instead, these two values h ere are adopted 
from the existing literature (IStairs et al.lll999^ to be regarded 
as a “reference standard”, we carefully adjusted the values of 
a and C, around the standard until the multi-wavelength pulse 
profile were reproduced. Though the final value of a = 40° 
is slightl y larger than the conve ntional values of a ;< 35° pre¬ 
sented in iDu et £0(120101120131) . but it is still in the allowable 
range. 
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Gaussian emissivities with distinct parameters for different 
region are also used to simulate profiles for three wavebands, 
in particular for the radio band. It is likely that radio emis¬ 
sion is originated in a non-caustic emission mechanism. Ra¬ 
dio emission from some pulsars might be generated in a in¬ 
coherent style. For instance, the radio wave could be divided 
into ordinary mode and extraordinary mode. When they are 
propagating in a pulsar’s magnetos phere, due to the r efraction 
effect for the ordinary-mode wave (IWang et al.l2014^ . the out- 
coming two-mode waves are separated which leads to a inco¬ 
herent and non-caustic emissivities. 

We predict that other pulsars’ multi-wavelength pulse pro¬ 
files with complex structures (e.g. PSR J2124—3358 with 
four or more radio peaks) and peak separations can be well 
reproduced by the self-consistent annular gap model with 
shifted emission altitudes. Do all normal pulsars (includ¬ 
ing young and millisecond pulsars) share the same physical 
mechanism for radio emission? This is a long-standing ques¬ 
tion that needs a larger pulsar sample and further studies. 
Joint modeling both radio pulse profile and polarization angle 
curve with considerations of various effects discussed above 
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4. CONCLUSION 

The multi-wavelength (radio, X-ray and 7 -ray) pulse pro¬ 
files of PSR B1821—24 can be explained by the annular gap 
model, and the emission region of the corresponding band are 
derived. The pulsed emission of all the three waveband are 
generated from the same magnetic pole. The shifted emis¬ 
sion altitudes for each waveband in the single magnetic-pole, 
which leads to a patchy emission pattern, can account for the 
puzzling properties of pulse profiles for PSR B1821—24. 
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